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Abstract. 

One of the very small number of serious alternatives to the usual concept of an 
astrophysical black hole is the "gravastar" model developed by Mazur and Mottola; 
and a related phase-transition model due to Laughlin et al. We consider a generalized 
class of similar models that exhibit continuous pressure — without the presence of 
infinitcsimally thin shells. By considering the usual TOV equation for static solutions 
with negative central pressure, we find that gravastars cannot be perfect fluids - 
anisotropic pressures in the "crust" of a gravastar-like object are unavoidable. The 
anisotropic TOV equation can then be used to bound the pressure anisotropy. The 
transverse stresses that support a gravastar permit a higher compactness than is given 
by the Buchdahl-Bondi bound for perfect fluid stars. Finally we comment on the 
qualitative features of the equation of state that gravastar material must have if it is 
to do the desired job of preventing horizon formation. 
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1. Introduction 

Although the concept of a black hole is well-established and generally accepted in the 
relativity, astrophysics, and particle physics communities, one sometimes encounters a 
certain amount of scepticism regarding the physical reality of the mathematical solution, 
and wariness regarding the interpretation of observational data £Q. The simplicity of 
the Schwarzschild solution comes at the expense of a central singularity, and an event 
horizon at the Schwarzschild radius -RschwarzscMd = 2M. Critics of black holes, at least 
the rational ones, seek an alternative configuration of matter that concentrates as much 
energy density as possible within a radius of R > 2M while avoiding the formation of 
the singularity and the event horizon. 

One of the small number of serious challenges to the usual concept of black holes 
is the "gravastar" (gravitational vacuum star) model that was recently developed by 
Mazur and Mottola El H|- I n the gravastar picture, or the very closely related 
quantum phase transition picture developed by Laughlin et al. E|, the quantum 
vacuum undergoes a phase transition at or near -RschwarzscMd where the event horizon 
would have been expected to form. Somewhat related models, differing in motivation 
and technical details, can be traced back to Gliner [2] and Dymnikova jHj. 

In the Mazur-Mottola model, a suitable segment of de Sitter space (with an 
equation of state p = —p > 0) is chosen for the interior of the compact object while 
the outer region of the gravastar consists of a (relatively thin) finite-thickness shell of 
stiff matter (p = p) that is in turn surrounded by Schwarzschild vacuum (p — p — 0).|| 
Apart from these three explicitly mentioned layers, the Mazur-Mottola model requires 
two additional infinitesimally-thin shells with surface densities a±, and surface tensions 
$±, that compensate the discontinuities in the pressure profile and stabilize this 5-layer 
onion-like construction, effectively by introducing delta-function anisotropic pressures 
[3 El II]- Since infinitesimally thin shells are a mathematical abstraction, for physical 
reasons it is useful to minimize the use of thin shells, either by successfully reducing the 
system to a 3-layer onion with one thin shell surrounded by segments of Schwarzschild 
and de Sitter space [TT], or more boldly (as in this article) by attempting to replace the 
thin shells completely with a continuous layer of finite thickness, 

Indeed, in the present article we demonstrate that the pressure anisotropy implicit 
in the Mazur-Mottola infinitesimally thin shell is not an accident, but instead a necessity 
for all gravastar- like objects. That is, attempting to build a gravastar completely out of 
perfect fluid will always fail. (Either the gravastar swells up to infinite size, or a horizon 
will form despite one's best efforts, or worse a naked singularity will manifest itself.) We 
derive this result by working with configurations where pressure is assumed continuous 
and differentiable, and analyzing the resulting static geometry using first the isotropic 

|| Because of the negative pressure core the gravastar exhibits some of the properties of the "tension 
stars" considered by Katz and Lynden-Bell [5], and by Comer and Katz ^Uj, though the intent and 
scope of the model is rather different. 
See also reference f° r related ideas. 
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TOV equation, and then the anisotropic TOV equation. 



2. The geometry 

We adopt coordinates that allow us to write any static spherically symmetric geometry 
in the form 

dr 2 



ds = — exp 



g(f) dr 



dt 2 + , w + r 2 dfi 2 . (1) 

1 - 2m{r)/r V 1 



We chose the two metric functions so that g(r) represents the locally measured 
gravitational acceleration, which is pointing inwards for positive g(r), and so that 
m(r) is the total mass-energy confined in a sphere with radius r. This interpretation 
can be justified by invoking the Einstein equations for the static stress-energy tensor 
= diag[— p,p r ,pt,pt\- The it-field equation then yields 

m(r) = An I p{r)f 2 df (2) 
Jo 

and the rr-field equation gives 

m(r) + Anp r (r) r 3 Aixr p(r) + 3p r (r) 
~ r 2 [1 - 2m(r)/r] ~ ~3~ 1 - 2m(r)/r ^ ' 

where we made use of the average density p{r) = m(r)/(|7rr 3 ). 

From now on, in the interests of legibility, we discontinue indicating the explicit 
r-dependence of all relevant functions. The remaining field equation for the transverse 
pressure is quite messy. Instead, we make use of the Bianchi identities and replace it 
with the covariant conservation of stress energy: 

- r = -{p + Pr)g+ • (4) 

dr r 

In the case of isotropic pressures p = p r = Pt this leads to the standard TOV equation, 

which can be written in any of the equivalent forms 

dp _ (p + p\g— (P + P) ( m + 4vrpr 3 ) _ An r (p + p) (p + 3p) 
dr r 2 [1 — 2m(r)/r] 3 1 — 2m(r)/r 

For anisotropic pressures we find it convenient to define the dimensionless anisotropy 

parameter 

A = (6) 
P 



so that in terms of this parameter the anisotropic TOV becomes 
dp r Airr (p + p r ) (p + 3p r ) 2 p A 



(7) 



dr 3 1 — 2m(r)/r r 

Before going on to explore some of the key features and consequences of equations 
© and (jZJ), we must define the class of spacetime geometries that we are particularly 
interested in. 
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3. Key features of smooth gravastar models 

To have a useful model, we should retain as much of standard physics as possible, while 
introducing a minimum of "new physics". In the spirit of Mazur and Mottola [21 EH E], 
and Laughlin et al. [HllEj, we will keep the density positive throughout the configuration 
but permit the pressure to become negative in the gravastar interior. To avoid 
infinitesimally thin shells one must then demand that the radial pressure p r is continuous 
(though the density need not be continuous, and typically is not continuous at the surface 
of the gravastar). Qualitatively the radial pressure is taken to be that of figure ^ 



\ J 














— ^ 

r 











r = r r max R 

Figure 1. Qualitative sketch of radial pressure as a function of r for a gravastar . 

That is, to "smooth out" the infinitesimally thin shells of the Mazur-Mottola 
gravastar model, we shall consider static spherically symmetric geometries such that: 

• Inside the gravastar, r < R, the density is everywhere positive and finite. 

• The central pressure is negative, p c < 0, and in fact p c = —p c . 
(We do not demand p = —p r = —p t except at the centre.) 

• The spacetime is assumed to not possess an event horizon. 
This implies that Vr we have 2m(r) < r. 

These three features, positive density, negative central pressure, and the absence of 
horizons, are the three most important features characterizing a gravastar. Other 
important features are: 

• To keep the centre of the spacetime regular, we enforce both p' r (0) = and 
Pc = Pr(0) = p t (0). 

• There should be a pressure maximum in the general vicinity of the Schwarzschild 
radius, r max ps -Rschwarzschiid, satisfying p r (r max ) > 0, and p' r (r max ) = 0. 

(This permits the physics in the region r ^> r max to be more or less standard.) 
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• There should be exactly two radii where the radial pressure vanishes: 

- The first pressure zero p r {r ) = 0, where p' r (r ) > 0, and 

— the second pressure zero p r (R) = 0, where p' r (R) < 0. 

The point R, (which by construction must satisfy R > r max > r ), is called the 
surface of the gravastar. 

• The pressure profile p r {r) should be continuous. 

(In contrast, it is sometimes useful to allow pt{r) to be discontinuous.) 

• The strong energy condition [SEC; p + p r + 2p t > 0] is definitely violated, at least 
near the centre of the gravastar. 

• We choose to enforce the null energy condition [NEC; p + Pi > 0] throughout the 
gravastar. In view of our first comment that density is everywhere positive, this 
implies that we are enforcing the weak energy condition [WEC; p + pi > and 
p>0]. 

• We impose no restriction regarding the dominant energy condition [DEC; p > and 
\Pi\ < p\i an d in fact we shall see that the DEC must fail in parts of the gravastar 
that are sufficiently "close" to forming a horizon. 

We comment that the Mazur-Mottola model might then be viewed as the limiting 
case where r max — > r , while p = —p r = —p t is strictly enforced for r < r , and where an 
additional thin shell is placed at the surface [so that p r (R~) > 0]. The 3- layer variant 
considered in is effectively the Mazur-Mottola model subject to the additional limit 
R — > r max, so in this model r = r max = R and there is a single thin shell at r max with 
de Sitter geometry inside and Schwarzschild geometry outside. 

Somewhat similarly, the Laughlin et al model may be viewed as the singular double 
limit r -> -Rschwarzschiid horn below, while R -> -Rschwarzschiid from above, and with 
p = —p r = —pt strictly enforced for r < -Rschwarzschiid- There is implicitly a singular 
infinitesimally thin shell located exactly at -Rschwarzschiid with infinite surface tension. 

The Gliner [2j and Dymnikova [Hj proposals all satisfy the constraint p = —p r 
everywhere throughout the configuration. For the purposes of this article we view this 
as an unnecessary specialization. 

One might wonder why we do not place the gravastar surface at r ? After all 
the pressure is by definition zero there so we can smoothly join it on to an exterior 
Schwarzschild solution. The reason for not doing so is a purely pragmatic one based 
on the fate of infalling positive-pressure matter. If we start without a positive-pressure 
region of type (tq, R), then any infalling positive-pressure matter that accumulates above 
r will automatically generate a positive-pressure region of type (r ,-R). The only way 
to avoid a positive-pressure region of type (ro, R) is if the negative-pressure matter 
in the (0, r ) region immediately catalyzes any infalling positive-pressure matter into 
negative-pressure matter. This scenario has its own risks, and for now we will keep the 
negative-pressure matter deep in the core, discretely hidden behind a layer of positive- 
pressure matter. (This is exactly what Mazur and Mottola did with their layer of stiff 
matter.) 
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4. Failure of isotropic pressure in the gravastar "crust" 

For the time being, we are only considering perfect fluids, so that the pressure is 
isotropic: p = p r = Pt and A = throughout this section. Consider first the isotropic 
TOV equation (JHJ) at r = r where the pressure is first zero. We find that 

dp 
dr 



4vrr pp (g) 



=ro 3 l-2m/r 

But the LHS is by assumption positive, while the RHS is by assumption negative. 
Therefore the isotropic TOV cannot hold at the point r . 

Secondly, consider the point of maximal positive pressure, r = r max . It follows from 
the isotropic TOV equation (J3J) that 

dp 



dr 



= 4?rr max (p + p)(p + 3p) 
3 l"2m/w ' 

But p > everywhere inside the gravastar, and by assumption p(r max ) > in the class 
of models we consider. So the LHS is zero while the RHS is negative. Therefore the 
isotropic TOV cannot hold at the point r max . 

In fact the objects under investigation have an increasing (radial) pressure for the 
entire range ro < r < r max while the (radial) pressure in the same interval is positive. 
That is the LHS of the isotropic TOV is positive in this region, while the RHS of the 
isotropic TOV is negative. Hence, it follows by contradiction from JSJ) that isotropic 
pressure in that interval is not able to satisfy the TOV equation, and thus a static 
spacetime geometry can only be obtained with anisotropic pressures. 

The same argument holds for a larger interval that extends below ro and into the 
negative pressure region: Assuming the NEC, we have p+p > 0. Therefore, by equation 
(JHJ) , we must have p + 3p < if the pressure gradient for isotropic pressures is to be 
positive. Now at the centre of the gravastar (p + 3p) c = p c + 3p c = —2p c < 0, while at 
r we have (p + 3p)o = po > 0. Therefore p + 3p changes sign somewhere in the interval 
r G (0, ro). Define r g to be the location where p + 3p changes sign. We conclude that 
pressure isotropy fails for the entire region where both p + 3p > and p' > 0, and quite 
possibly fails for an even larger region. That is, isotropy definitely fails on the region 
r g < r < r max , a region which definitely encompasses r < r < r max . 

Let us call the region r e (0, r g ) the "core" where the physics is qualitatively (if 
not necessarily quantitatively) similar to that of de Sitter space. In particular in the 
"core" the local acceleration due to gravity [given by equation ©] is outward. Similarly 
let us call r e (r 9 ,r max ) the "crust", where physics is still definitely "unusual". In the 
"crust" the local acceleration due to gravity is inward, but the pressure still rises as 
one moves outward. Finally, call r G (r max , R) the "atmosphere", where the physics is 
"normal" , or rather as normal as it is going to get in a gravastar. In the "atmosphere" 
the local acceleration due to gravity is inwards, and the pressure decreases as one moves 
outwards. (See figure El) With these definitions, we see that pressure is guaranteed 
to be anisotropic throughout the "crust". Note that even if we were to dispense with 
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the entire positive-pressure region by chopping the gravastar off at r , there is still an 
anisotropic crust in the region (r g ,r ]. 



Pr A- 



core 



crust atmosphere 




Figure 2. Qualitative sketch of gravastar labelling the "core" 
"atmosphere" . 



"crust" , and 



To conclude this section, we summarize that a static spherically symmetric object 
with positive density, negative central pressure and vanishing pressure at the surface 
cannot be supported by isotropic pressures alone — there are no perfect fluid gravastars. 



5. Schwarzschild interior solution and the Buchdahl— Bondi bound 

A particular illustration of what goes wrong when one tries to build a perfect fluid 
gravastar is provided by the Schwarzschild interior solution. In this model one assumes 
constant positive density p* throughout the whole star and invokes the isotropic 
TOV equation Despite common misconception, this does not mean the star is 
incompressible or that the speed of sound goes to infinity. See for instance [T31 p. 609 
ff] and JT3]. In this situation one has the well-known analytical solution [T3] 

y/1 - 2m*(r)/r - ^1 - 2M/R 



^ P * 3y/l - 2M/R 2m,(r)/r 

where M and R are the total mass and surface radius and m*(r) 
central pressure 

1 



(10) 

M(r/Rf. The 



Pc = P*: 



2M/R 



'3y/l - 2M/R- 1 

certainly diverges for 2M/R ^8/9 and is positive for all 2M/R < 8/9. This is a first 
indication of the existence of the Buchdahl-Bondi bound, but by itself is not enough to 
derive this bound for arbitrary perfect fluid spheres ^H] • Since we are only interested in 
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this model as an example, we need not worry about fully general statements and simply 
note that 

8 2M 

Pi < « 



9 < 



R 



< 1. 



(12) 



While this is a perfectly sensible solution in the mathematical sense one would generally 
rule it out physically because of the negative central pressure. If we leave our prejudice 
against negative pressures aside, adopting the gravastar philosophy, we find that with 
the choice of ()12|). the pressure profile ()10j) will have a first order pole at 



^pole 



3R\ 



1/9 



< R. 



(13) 



2M/R 

Note that as 2M/R goes from 8/9 to 1, the position of this pole moves from the centre 
of the star to the surface of the star. The situation is qualitatively sketched in figure H3 
Note that the NEC must be violated sufficiently close to the pole. 



Pr A 



Pc 




r = 



Figure 3. Qualitative sketch of pressure pole in the interior Schwarzschild solution 
for 2M/R > 8/9. 

This is now clearly unphysical, even if one is willing to accept negative pressures 
and even violations of the NEC. It is unphysical for more prosaic reasons because the 
pressure pole implies a curvature singularity — in fact the Rmt an d Rieie orthonormal 
components of the Riemann tensor are infinite, so that one has a naked singularity. 
The reason why we mention this specific example is because we shall soon see this 
"pressure pole" behaviour is generic — continuous solutions with isotropic pressure are 
not possible as we have already shown in section |3J + 

+ Consider for example reference where a perfect fluid Chaplygin gas p oc 1/p is considered. The 
surface of their configuration occurs at p = where p = — oo, at least as one approaches the surface from 
below. So the surface of their configuration is a naked singularity, in agreement with the observations 
above. 
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6. The fate of a negative pressure perfect fluid sphere 

Consider now an arbitrary gravastar (of the type defined above). We note that the 
"compactness" x — 2m(r)/r satisfies 



2m(r) 



2m 8n r . 8ir r . . , . . . 

l7r p r ~ -r = — ftp ~ p\ = — [ 3 (/° + p) - (p + W ( 14 ) 



But the first term on the RHS is non-negative by the NEC, while the second term is 
by definition negative on [0, r g ), so the compactness 2m(r)/r is monotone increasing on 
the range [0,r g ). This is more unusual than one might expect from the simplicity of 
the argument. While the compactness of a normal perfect fluid star tends to increase 
as one moves outwards, in a normal star it also can be subject to oscillations that make 
the overall picture quite subtle [13 HHj . 

Consider now a perfect fluid sphere with negative central pressure that satisfies the 
NEC. Since we have already seen that isotropy is violated on (r g , r max ), the only way we 
can maintain the perfect fluid nature of the sphere is if r g — > oo. (Whence also r — > oo.) 
But since we do not want a horizon to form, the compactness must be bounded above 
by unity. And since we have just shown the compactness is monotonic we now see 

Um 2m(r) 6(0,1]. (15) 

r-^oo T 

So not only does a NEC-satisfying perfect-fluid gravastar expand to infinite volume, it 
also has infinite mass. 

To avoid the physical size of the gravastar blowing up to r — > oo, our options now 
are rather limited: We could permit the development of a horizon at finite r, which 
defeats the whole point of the exercise, or we could permit something even worse. If we 
permit NEC violations then it is possible to arrange for the development of a pressure 
pole at finite r < r . To see how this is possible, consider the configuration 

r 

p(r) 7n , (16) 

r — r p 

and let us check that it is compatible with the isotropic TOV equation. First for r < r p 
the pressure is negative, and for r > r p the pressure is positive, which is appropriate for 
a modified gravastar model. Secondly 

^(^<°- < 17 > 

Third, assuming p remains finite, close to the pole (r w r p ) we have 

(p +p)(m + Anpr 3 ) ^ [T / (r - r p )} [4irr 3 T / (r - r p )} 

r 2 (l — 2m/r) r 2 (l — 2m{r)/r) 

Airr V 2 1 

v < 0. (19) 



(1 - 2m(r p )/r p ) (r - r p ) 2 
So the TOV equation can be satisfied in the vicinity of the pole provided we set 
l-2m(r p )/r p 



4nr p 



> 0. (20) 
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That is, if the gravastar configuration is perfect fluid and finite in extent — then this 
pole is the only way the TOV equation can be satisfied. (It is easy to check that higher- 
order poles do not even have this nice property of being compatible with the isotropic 
TOV equation.) Now despite the fact that the presence of a simple pole in the pressure 
is compatible with the isotropic TOV equation, we must reiterate that such a pressure 
pole is unphysical because it is a naked singularity — the physically correct deduction 
form this analysis is that gr avast ar-like objects must violate pressure isotropy. 



7. Bounds on the pressure anisotropy 

Once we accept that perfect fluid spheres are not what we are looking for to model 
gravastars, one might wonder what happens to the Buchdahl-Bondi bound for isotropic 
fluid spheres. It has been shown that for p' < and p t < p r the 8/9 bound still holds. 
However if the transverse stress is allowed to exceed the radial stress, pt > p r , then the 
upper limit shifts to 2M/R < k < 1, where k depends on the magnitude of the maximal 
stress anisotropy ^Hj-* In the gravastar picture, we shall soon see that in the crust 
p t > p r , and that the compactness of a gravastar is not limited by the Buchdahl-Bondi 
bound, but only by the magnitude of the maximal pressure anisotropy and the regularity 
of the metric, i.e. 2m(r)/r < 1. Let us now make these qualitative considerations more 
quantitative. To do this, let us rewrite (J7J) in the form 

As *z» 'r* + ( 1 + *) 9 l. (21) 

P 2 LP \ P J . 
In section |3] we have determined the smallest interval in radii for which anisotropic 
pressure is necessary to be (r g , r max ]. We can now ask for explicit bounds on A in that 
interval. By inserting the definition of g, (|7j). into (|21|). we get 



A = ^ 



Pr ( p r \ m + 4irp r r 3 



1 + £_l ) _ 1 ~" r ' " (22) 

2 |_ P ' V p / r 2 [1 - 2m/r] 

For the interval r G [ro, r max ], by making use of p' r > and p r > 0, we find the simple 
lower bound 

1 2m lr , , 

A > - > 0. 23 

~ 4 1 - 2m/ r y ' 

Now in the region [r ,r max ] we have p t > p r > 0, consequently if the DEC is to be 
satisfied we must at the very least have A < 1. But this is guaranteed to be violated 
whenever 2m /r > 4/5. 

That is: If the gravastar is sufficiently close to forming a horizon, in the sense that 
2m/r > 4/5 somewhere in the range [ro,r max ], then the DEC must also be violated at 
this point. Jj Consequently, any gravastar that is sufficiently close to forming a horizon 
will violate the DEC in its "crust" . 

* For related work on anisotropic stars, see for instance |20| . 

jj And even if we were to discard the entire positive-pressure region (ro,i?), we can nevertheless still 
apply this bound at ro itself: If 2m(r^)/ro > 4/5 then the DEC is violated at Vq. 
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For the interval r £ {r g ,r ) we find the considerably weaker bounds 

rp' . rp' 1 2m /r , . 

- 2p 2p 4 1 - 2m/r v ; 

where we have used p' r >0,p r < 0, and the NEC. 
8. Minimizing the anisotropic region 

Let us now attempt to minimize the region over which anisotropy is present. It is easy 
to see that at r g we have 



A(r g ) = ^ > . (25) 

Therefore, in the case where the anisotropy is confined to the smallest interval possible, 
we want p' r {r g ) = 0, corresponding to an inflexion point for the radial pressure. At the 
point r of zero radial pressure, the anisotropy cannot vanish: 

A(r ) = ^ + 1 2m/r ° >Q. (26) 
v u; 2p 4 1 - 2m/r V ; 

At the point of maximal radial pressure, the anisotropy also has to be non-zero, at least 
if we take the limit from below: 

A(r„J = i 2m '[T fl + ^Ul + ^M>0. (27) 



4 1 - 2m/r max \ p J \ m 

Beyond the peak Vr > r max , it is possible to arrange A = 0, though at a price: If 
we wish to confine the anisotropy to the smallest interval possible, we have to set 
A(r — > r+ ax ) = which leads to a discontinuity in p' r and A as well as a "kink" in the 
pressure profile p r at r max . (However, p r itself is still continuous, as is the density p.) 
Indeed we then have 

p' r (rL x ) = — A(r max ). (28) 
r„ 



max 



The implications of confining the pressure anisotropy to the smallest interval possible 
are shown in figure IH where the anisotropy is confined to the region r £ (r 9 ,r max ]. 

The Mazur-Mottola model is now recovered as the limiting case where r g — > 
''max, and so all the important anisotropy is confined in their inner thin shell. 
The anisotropy A — ► oo, because p' — > oo. Effectively A is replaced by choosing an 
appropriate finite (in this case negative) surface tension & and surface energy density a 
which is given by the Israel-Lanczos-Sen junction conditions [TT| I21j. 

The second, outer thin shell which is present in the Mazur-Mottola model is 
not a physical necessity, but is merely a convenient way to avoid an infinitely diffuse 
atmosphere that would arise otherwise from the equation of state p = p. A finite surface 
radius R can be modelled by altering the equation of state slightly to include a finite 
surface density p$ which is reached for vanishing pressure: p(p) = ps + P- Then, the 
outer thin shell can be omitted when joining the gravastar metric onto the Schwarzschild 
exterior metric. 




Figure 4. Qualitative sketch of radial pressure as a function of r for a gravastar that 
is as near as possible to a perfect fluid. Note the inflexion point at r g and the kink at 

^max • 



9. Features of the anisotropic equation of state 

In the case of a perfect fluid, the geometry is completely defined by the set of differential 
equations (j2j), (jHJ) plus initial conditions and an equation of state, which is commonly 
written as p = p(p) or equivalently p = p(p). 

Once we accept the need to abandon the notion of a perfect fluid, we have 
to replace the isotropic TOV equation (JHJ) with its anisotropic counterpart (jJJ) and 
thereby introduce an additional (free) function A(r). To close the set of equations, it 
is now necessary to define two equations of state. Naively, one might choose p r {p) 
and Pt{p), which is a rather strong assumption that forces p r and p t to change in 
lock-step. Alternatively, one could simply postulate a density profile p(r) (which is 
the strategy adopted, for example, in [Hj) or equivalently a pressure profile, plus one 
equation of state. One could also (rather unphysically) choose to specify any two 
profiles p(r), p r (r), and Pt(r) by hand, and use the anisotropic TOV to calculate the 
remaining profile. Yet another possibility to obtain a well defined solution lies in finding 
an additional (differential) equation, that might be motivated by some appropriate 
variational principle, for example by minimizing the "total anisotropy" 



1 



R 



R 



(pA)Mr 5 



R 



0. (29) 



(pA) 2 dr 

No matter whether it is an equation of state for the transverse pressure or a variational 
principle for the total anisotropy, the extra equation should (if the gravastar model is to 
be even qualitatively correct) be responsible for maintaining stability of the gravastar 
over a wide range of total masses and central pressures. In other words, the closed set 
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of equations must be self-regulating if it is to be physically interesting — the gravastar 
should shift to a new stable configuration when the total mass changes. 

The central point of this article is that if you wish to believe in gravastars you must 
accept that there will be regions of anisotropic pressure where A > 0. This suggests 
that it might be most efficient to choose the two distinct "equations of state" as being 
equations for p and for A. But what variables should these equations of state depend 
on? An obvious candidate is the radial pressure, but in view of the inequality (|23|) it is 
clear that the gravastar material, if it is to succeed in generically avoiding the formation 
of horizons should be sensitive to the "compactness" 2m(r)/r. That is, we should posit 
equations of state of the form 

p = p(r,p r ,2m/r); A = A(r,p r , 2m/r); (30) 

and (making everything as explicit as possible) solve the paired differential equations 

dp r (r) [p(r,p r (r),2m(r)/r) + p r (r)] [2m(r)/r + 8n p r (r) r 2 ] 



dr 2r [1 — 2m{r)/r 

dm(r) 



t 2 p(r,p r (r),2m(r)/r) A(r,p r (r),2m(r)/r) 



4n p(r,p r (r),2m(r)/r) r 2 . (32) 
dr 

Whether or not the gravastar model ultimately succeeds in its goals depends on whether 
or not one can find physically realistic equations of state that have the effect of bounding 
2m(r)/r < 1 for large ranges of p c and total mass. 

Now there is an important issue of principle here: Traditional relativists are 
somewhat nervous when the variable 2m/ r enters the equation of state, arguing that 
2m I V is not detectable by local physics, and that an equation of state that depends on 
2m jr somehow violates the Einstein equivalence principle. This is not correct for the 
following reason: 2m jr is certainly measurable by quasi-local effects in small but finite 
size regions. To see this note that for any static spherically symmetric spacetime the 
orthonormal components of the Riemann tensor are (see, for instance, [221 P HO]): 

4tt 

RtHf = y [3(P " Pr + 2p t ) - 2p\ (33) 

4:71 

Rme = Y [3pr+ ^ (34) 
4ir 

Rfm = Y [3p ~ p]] (35) 
4-k 

Reiei = [2p]- (36) 



3 

The point here is that the Riemann tensor is certainly measurable in finite-sized regions, 
so in particular p is measurable. Likewise r is measurable in finite-sized regions, and 
therefore 2m jr = (87r/3) p r 2 is measurable. Consequently the compactness 2m jr is 
a quasi-local measurable quantity and it can meaningfully be put into the equation 
of state without violating the equivalence principle. (In either the Mazur-Mottola 
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scenario (21 El Hj or the Laughlin et al. scenario EH the gravastar material is assumed 
to be a quantum condensate, and therefore sensitive to non-local physics. The point of 
the current discussion is that we do not need to appeal to quantum non-locality to get 
2m/ r into the equation of state — properly understanding the equivalence principle is 
enough.) Note that while this argument demonstrates that possible equations of state 
are at least conceivable, that is not the same as explicitly demonstrating that such 
equations of state actually exist. That is a challenge which we leave for the future. 

10. Results and Discussion 

In this article we have delineated the qualitative features one would expect from a 
gravastar configuration that avoids delta function transition layers and has finite regular 
pressure profile at all locations. We have used these qualitative features to place 
constraints on the anisotropy parameter A, demonstrating that perfect fluid gravastars 
are a lost cause, and extracting some generic information regarding what the gravastar 
equation of state should be. 

Specifically, finite-sized gravastar- like objects (with the key defining feature being 
negative central pressure) must exhibit anisotropic stresses in their "crust", which is 
the region where the pressure is increasing as one moves outwards but the local force of 
gravity is inwards. Trying to build a perfect fluid gravastar results either in an infinite- 
size infinite-mass object, or in a naked singularity as the pressure exhibits a simple 
pole. Assuming the WEC, the magnitude of anisotropy required in the crust can be 
explicitly bounded in terms of the local compactness 2m(r)/r, and becomes arbitrarily 
large for gravastars that are sufficiently close to forming a horizon, in which case the 
DEC must be violated. Consequently, if one demands that the equation of state for 
gravastar matter results in configurations that are horizon-avoiding for large ranges of 
total mass and central pressure, then it is difficult to avoid the conclusion that the 
equation of state must depend on the local compactness 2m(r)/r. This is certainly an 
unusual equation of state, but we emphasise that (when properly understood) this does 
not violate the equivalence principle. 

While we are personally agnostic as to the existence or non-existence of gravastars, 
we feel it is important to understand what their general properties might be in order to 
have a clear understanding of what the observational evidence regarding astrophysical 
black holes is actually telling us. 
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